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Abstract 

The effects of Rashba spin-orbit coupling on BCS-BEC crossover, the 
condensate fraction and pair coherence lengths for a two-component 
attractive Fermi gas in two dimension are studied. The results at T = 
OK indicate that (1) when the strength of SOC is beyond a critical value, 
BCS-BEC crossover does not happen in a conventional sense; (2) SOC 
enhances the condensate fraction, but suppresses pair coherence lengths. 



PACS number(s):03.75.Ss, 05.30.Fk, 67.85.Lm 

In a crystalline solid spin-orbit coupling (SOC), which occurs naturally in systems with 
broken inversion symmetry and makes the spin degree of freedom respond to its orbital mo- 
tion, is responsible for many interesting phenomena, such as magnetoelectric effect [ll[2l[3], 
visionary Datta-Das spin transistor [HIS], topological insulator [HIIT] and superconductivity 
[8l[9]. Taking topological superconductivity for example, it has been predicted to occur in 
superconductors with a sizable spin-orbit coupling in the presence of an external magnetic 
field [ini im HI [13l HI] . In these systems the transition to topological phases requires that 
critical magnetic field is much larger than the superconductivity gap above which an s-wave 
superconductor is expected to vanish in the absence of SOC. It is SOC that competes with 
a strong magnetic field to give rise to a topological superconducting phase. 
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It is widely known that ultracold atom systems can be used to simulate many other 
systems owing to their many controllable advantages and operabilities [15\ [161 Cer- 
tainly the simulations to SOC, which are generally equivalent to produce non-abelian gauge 
potential with optical |18l 1191 1^0] or radio- frequency fields |21j . are also possible and have 
been realized in a neutral atomic Bose-Einstein condensate (BEC) by dressing two atomic 
spin states with a pair of lasers [22]. Motivated by such a pioneer experiment and a practical 
proposal for generating SOC in '^^K atoms [23], BCS-BEC crossover in the two-component 
Fermi gases with SOC have been widely studied [SI 113 [261 [271 [281 IMl [3Ql EI] ■ On the 
one hand for balanced case SOC not only leads to an anisotropic superfluid [24], the sig- 
natures of which could be observed in the momentum distribution or the single-particle 
spectral function of atomic cloud, but also significantly enhances the superfluid transition 
temperature when scattering length < 0, while suppresses it slightly when > [25J. 
In addition by adjusting the strength of SOC, one can engineer a BCS-BEC crossover even 
with a very weak attractive interaction that is unable to produce a two-body bound-state in 
free vacuum [26]. On the other hand for imbalanced case SOC and population imbalance are 
counteracting, and this competition tends to stabilize the uniform superfluid phase against 
the phase separation. However, SOC stabilizes (destabilizes) the uniform superfluid phase 
against the normal phase for low (high) population imbalances [30j . 

In this paper we consider BCS-BEC crossover of two-component Fermi gases with SOC 
positioned in a two-dimensional (2D) space, and are interested in the evolutions of con- 
densate fraction and pair coherence lengths along the crossover. As is known to all that, 
without SOC, a spin-up fermion pairs with a spin-down fermion, i.e. fermion pairs happen 
in the singlet channel. This leads to that the condensate fraction comes from the contri- 
bution of singlet pairs [32l [33] and the coherence length of singlet pairs is defined [3ll |35] • 
In the presence of SOC, SOC also induces triplet pairs in the system besides singlet pairs. 
Thus at this time both singlet and triplet pairs contribute to the condensate fraction, and 
we must define two coherence lengths related to singlet and triplet pairs. 

The Hamiltonian of the system we consider is 



H= /dV|^M/t^(r 



-Q 

— + VR{axPy - ayPx) - /i 



*^(r) - U^{f)^Uf)^'^{f^^Ur) \ , (1) 



where a Fermi atom of mass m for spin a is described by the field operator ^ci(r). ax 
and ay denote the Pauli matrices in the x and y directions. —U{U > 0) corresponds to 
attractive contact interaction among fermions and fi is the chemical potential. Without loss 
of generality we assume SOC to be Rashba type and its strength to be vr. Transforming 
the field operator ^^(r) into momentum space 

^a{f^ = ^j:'^a{k)e^'-', (2) 



where V is the volume of the system, and introducing the mean-field order parameter 
^ = yY.k < *4.(-^)*t(^) > Hamiltonian ^ is written into 

H = 5]efe[^{(^)*tW + ^l(^)^4W] + ^'fc^t(^)^4(^) + ^'fe^Iw*t(^) 



-A*{(^)^|(-^) - A^^{-k)^^{k) + (3) 
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with ^fc 



2m 



fj, and Qk = fiVR{ky + i/c^ 



Explicitly the Hamiltonian ([3]) is second order about field operators and can be solved 
exactly. To attack this goal we choose to use imaginary time Green function method [36j 
since in this frame some interesting physical quantities, such as atom number and order 
parameter, can be directly deduced from Green functions. Defining two normal Green 
functions G^tW = - < Tr'i'-i{k,T)^{k,0) >, G^,^{t) = - < r^^'|(fc, r)^'|(A;, 0) > and 

two anomalous Green functions F|^(r) = 



< Tr^li-k,T)^ik,0) >, F^^ir) = - < 
Tr'^^^{—k,T)'^^{k,0) >, and using the time evolution of an imaginary time operator O, 
= [0,H], we obtain the equation of motion 



M4x4(Gtt, nt' Pnf = (1' 0, 0, 0)^ 



(4) 



where G^^, Gj^-f-, F^-^, F-^^ are Fourier transformation of G^^{t), G|-|-(r), -F4.|(r), F|-f-(r) respec- 
tively and 
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(2ra+i)7r pgpj.gggjj^^jjg ^]2g Matsubara frequency. Prom ([5|) the quasiparticle ex- 



with Wn 

citation spectrum is obtained by replacing iwn with w and setting detM4x4 



I3h 

0, leading 

to i?^ -p = {^k T \Qk[f' + Naturally for vr = 0, the excitation spectrum reduces to the 
standard BCS form = ^? + A^. By using E^ zr. Green functions can also be figured out 
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From ([6]) we know for a general k, Green function F^^ is non-zero. This signifies that SOC 
induces triplet fermion pairs although the interaction is pure s-wave. 

The evolution of the system is completely decided by the equations of order parameter 
A and particle number 



- = -E 
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(7) 
(8) 



by using the relations A = 'Ek,iw„ -^it and iV = 'Ek,iw„ Gn- I* is easily found that 
the equation ([7]) is divergent due to the fact that we use contact interaction to simulate the 
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true two-body potential, whose Fourier transformation should fall off at large momentum. 
To regulate the divergence the strength of contact interaction U should be substituted into 



l = ly I (9) 

where is the 2D two-body binding energy |37j . 

We self-consistently solved the equations d?]) and ^ for the different strengths vr of 
SOC at T = OK. In Fig. 1(a), the evolution of the order parameter A is shown and we 
can find that the larger vji is, the larger A is. From this perspective the existence of 
SOC enhances the superfluidity of the system. By comparison the evolution of /x is more 
interesting and plotted in Fig. 1(b). It is easily found that with the increase of vji, the 
chemical potential n becomes negative even if the two-body binding energy is very 
small. Hence from a conventional viewpoint, that BCS-BEC crossover is exactly realized 
when the chemical potential /i crosses zero, this fact suggests that when vr is beyond a 
critical value Vj^, BCS-BEC crossover does not happen at all and the system is into BEC. 
From numerical work vj^ w vp/\/2, where vp is 2D Fermi velocity and has a relation with 
particle density n vp = hyj2'Kn/m. It is to be noted that 2D Fermi energy ep = n-nfi^ jm 
and Fermi wavevector hp = \/2vrn. 

Now we analyze the condensate fraction of fermions in the BCS-BEC crossover with 
SOC. In terms of a superfluid Fermi system the condensate fraction no is related to off- 
diagonal long-range order (ODLRO) [38] and corresponds to the maximal eigenvalue A^O; 
divided by the whole particle number iV, of two-particle density matrix 

/02(ricJi,f2CJ2 : ficJl,f^cT^) =< ^'t^(fi)^'t^(f2)^',. (fa)^'^. (fl) > . (10) 

According to Leggett's book [39j, A^'o can be decided as follows 

iVo= E / / dr2|$(riai,r2a2)|2, (11) 

0-1,0-2 

where ^(ncji, r2<72) =< ^ axif"^^ aiif2) > is an anomalous average which arises as a result 
of spontaneous breaking of the U{1) gauge symmetry. 

According to the above theory and a fact that SOC induces triplet pairs, the condensate 
fraction in the presence of SOC is 
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J2 [\ < ^t(^)^;(-^) > I' + I < *t(^)*t(-^) > I'] • (12) 



The pre-factor 2 comes from the contributions of spin summation and time-reversal symme- 
try. In contrast to the case without SOC, there is an extra contribution to the condensate 
fraction from triplet pairs. Following the same procedure deriving the equations ([7]) and 
([8]), we have 
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When vn = 0, -Bfc,- = Ek^+ = Ek and uq = ^X^fctanh^ ^Ek/E^, same as the results in 
[33j . In Fig. 2 we calculate the condensate fraction no for different vji at T = OK. The 
results are twofold. Firstly SOC also enhances the condensate fraction, which is consistent 
with Fig. 1(a). In addition, with the increase of rapidly increases to a large value for 

a large vpt- Maybe this phenomenon can be illustrated from the results of Fig. 1(b), that 
for a large vr there is not BCS-BEC crossover and the system is situated in BEC, building 
on the fact that for a weakly interacting BEC, almost all atoms are into the condensate. 
A recent paper [lOj also calculate the condensate fraction in the same system but do not 
think over the contribution from triplet pairs, so their result is qualitatively incorrect. 

At last we determine the coherence lengths for singlet and triplet pairs. On general 
ground, information on pair coherence lengths can be extracted from the pair-distribution 
function 

9.^,a,{r) = ^< ^t^(r)^'t^(0)M/.,(0)^.,(r) > . (14) 

Following the same spirit that the Hartree-Fock term has been neglected in the Hamiltonian 
([3]), at the mean-field level ([H]) becomes 



g.^,a,ir) = ^\< ^Uf^^UiO) > |2, (15) 



and pair coherence lengths Cnair' ^iair singlet and triplet pairs can be obtained as 



with (pf. — ^ 
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Without 



SOC, ifl = but ipl = ^tanh^Ek/Ek, consistent with the results in [Ml ES]. Fig.3 
and Fig.4 describe the behaviors of (,pair and Cpair process of evolution at T = OK, 

respectively. Very explicitly, SOC suppresses pair coherence lengths for both singlet and 
triplet pairs, and for a large vr, pair coherence lengths rapidly decrease to a small value. 
We think that this phenomenon can also be understood from Fig. 1(b) in the light of the 
fact in BEC region pair coherence length is much smaller than in BCS region. Besides by 
comparing such two figures, triplet pair coherence length Cpair always is larger than singlet 
pair coherence length ^p^j^. Physically this is the result from Pauli exclusion principle. 

In summary we have discussed the evolution from BCS to BEC superfluids in the pres- 
ence of Rashba SOC in two dimension and shown that when the strength of SOC is beyond 
a critical value, BCS-BEC crossover does not happen in a conventional sense. In addition, 
we studied the evolutions of the condensate fraction and pair coherence lengths. The results 
indicate that SOC enhances the condensate fraction, but suppresses pair coherence lengths. 
Furthermore we also give some physical interpretation for some phenomena. 
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Figure 1: The evolutions of order parameter A (a) and chemical potential n (b) as a 
function of two-body binding energy es, in units of the Fermi energy ep, for different the 
strengths of SOC at T = OK. In (a) from bottom to top and in (b) from top to bottom 
vr/vf = 0,V2/2,1,V2,2. 
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Figure 2: The evolution of the condensate fraction no as a function of two-body binding 
energy es, in units of the Fermi energy ep, for different the strengths vr of SOC at T = OK. 
From bottom to top vr/vf = 0, \/2/2, 1, ^/2, 2. 
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Figure 3: The evolution of singlet pairs coherence length ^^^^^ as a function of two-body 
binding energy ef?, in units of the inverse Fermi wavevector kp and Fermi energy ep 
respectively, for different the strengths vr of SOC at T = OK. From top to bottom 
vr/vf = 0,V2/2,l,^/2,2. 
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Figure 4: The evolution of triplet pairs coherence length Cpaij- ^ function of two-body 
binding energy es, in units of the inverse Fermi wavevector kp and Fermi energy ep 
respectively, for different the strengths vr of SOC at T = OK. From top to bottom 
vr/vf = V2/2,1,V2,2. 
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